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Abstract 

Within the QCD factorization approach, we calculate the process- and polarization-dependent non- 
factorizable terms a a of the J/^ + K* decay. The longitudinal part oq is infrared convergent and 
large enough to agree with recent experimental data, provided that the B-K* form factors ^i(m^) 
and A2(m|,) satisfy some constraints met by some (but not all) models. The transverse parts a± on 
the other hand are infrared divergent, the procedure used to handle such divergence is discussed in 
relation with the B^ J/^' + K case in which the same problem arises. Our nonzero phases of the 
helicity amplitudes are consistent with experimental data recently measured for the first time by the 
CDF and BaBar groups. 
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1 Introduction 



Among a hundred hadronic decay modes of the B mesons[[I]] investigated from both experimental and 
theoretical sides, the process J/^ + K*(892) is particularly interesting for many reasons: 

(i) - It is the first color-suppressed B decay observed in 1994 by the Argus group[^] with the largest 
branching ratio for its class. Since then, important measurements are intensively explored in great 
detail by the Cleo[§], CDF[§], BaBar[|] and Belle[g] collaborations. 

Using both angular and time distributions to separate the CP-even from the CP-odd eigenstates, the 
asymmetry between B^J/^f + K* and B^J/^' + K - which directly gives sin 2/? at the same 
degree of precision as its companion "golden" mode B(B) J/^f + Ks - is central to our under- 
standing of the standard Kobayashi-Maskawa CP violation mechanism. Furthermore, new physics 
beyond the Standard Model can be hinted by consistently comparing the (3 angle obtained from this 
CP asymmetry with its unitarity triangle value determined by other experiments (the e measurement 
in K decays, B'^-B mixing, Vub/Vcb---)- Are they unambigously equal ? 

(ii) - Over both the vector -i- pseudoscalar B— > V -i- P and the pseudoscalar -i- pseudoscalar B^ P -i- 
P modes, the advantage of the vector -i- vector B^ V -i- V decays (B^ J/^* + K* considered here) 
stands out in the possibility of detailed analyses of the three helicity amplitudes. The three decay 
amplitudes (one longitudinal and two transverse) denoted by Hq, H_i and H+i could be separately 
determined both for their magnitudes |H;^| and phases 5\. These and 5\ analyses provide a pow- 
erful tool to test not only the naive factorization method[^ - usually adopted to deal with exclusive 
two-body hadronic decays - but also the real and imaginary part of the nonfactorizable terms which 
are calculable in QCD approaches [|8], ^ can then be confronted with experiments. We note that such 
analysis cannot be done in B^ V-i-P and B^ P-i-P decays since only the absolute value of a single 
amphtude (the equivalent of |Ho|) can be measured in these processes. 

(iii) - The transverse H„i and H+i amplitudes with both magnitudes and phases provide also an useful 
way to test robustness of factorization manifested through the V ^ A property of the effective weak 
currents in the Wilson operator product expansion (OPE). It would answer the question[[ro|] whether 
or not the interactions between the hadronic decay products, usually called final-state interactions 
(FSI), are strong enough to flip the quark spin in color-suppressed B decays. Although intuitively this 
spin flip unlikely occurs, this possibility could be tested however. 

(iv) - Improved by QCD which gives the ag corrections to the decay amplitudes in OPE, the three 
"heavy to Ught" B-K* form-factors in B^ J/^* + K* : Ai(g^), ^2{q^) and V(g^) can be determined 
and compared to models given in the literature. These form factors are useful for other decays, in 
particular B^ p -i- K* and B^ $ + K*. 

Motivated by new experimental data[Q, and recent theoretical developments [|8|, ^, we are trying in 
this paper to investigate some aspects of the B^ J/^' + K* process. 



2 Decay amplitudes 

2-1 Generality, Polarizations and Angular Distributions 

The most general B^ Vi + V2 helicity amplitude takes the following form in which we adopt the 
sign convention of[pl|] 

H,(b(P) ^ V(j.i,ei)+Vfe,e2)) =e;'^(A)er(A) ( g^.uA + + \e,,^p^^c\ , 
\ y ± / / 1 r- mim2 mim2 / 

(1) 

where A stands for the three helicities 0, ±1 of the massive vector mesons with polarizations e^(pi), 
€2(^2)- Since the initial B meson is spinless, the two final vector mesons share the same helicity A. 
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Here M, mi, 1712 are masses of the B, Vi, V2 mesons with four-momenta P,pi,p2 respectively. The 
A and B associated to the S and D waves are CP-even, while C corresponding to the P wave is CP- 
odd. We note that in two-body decays, the Lorentz invariant amplitudes H;^ have a mass dimension, 
so are the quantities A,B,C. From (1), we get 



Ro = -(aA + {a^ -l)Bj , B.±i = A ± Va^^ C , (2) 

with|0] 

pi ■p2 M"^ -mj- ml 
a = = . (3) 

m\m2 2mim2 

Also we define Kc as the common momentum of the outgoing mesons Vi (or V2) in the B rest frame: 

2 ^ A(M^mf,mi) 
4M2 

where X{x, y, z) = + + — 2{xy + xz + yz). Thus ^/a? — \ = MKc/ (mim2). 
From (1) and (2), the decay rate in each polarization state is given by: 

= , (4) 

with r = XIa is the full decay width r(B ^ J/^' K*). 

In the following we normalize the partial widths Tx = Tx/T for the three independent polarization 
states 

Ha 



such that Ea |Ha P = 1- 

The normalized dimensionless spin amplitudes Aq, A\\ and A_i_ are related to the helicity ones by 

/I - fr A - H+i +H-1 H+i - H_i 
^0 — -no , ^11 — , ^± — , 

with again \ Ax P = 1. To proceed to the Aq, A\\,Ax_ determinations, angular measurements are 
necessary. For that, let us define the transversity angles 6'tr and <I>tr as the polar and azimuthal angles 
of the £^ descended from J/^* ^ + £~ decay in the J/^' rest frame. The K* helicity angle Ok* 
is the angle between the K meson direction (coming from K* vr-i-K) and the opposite direction of 
the J/^' in the K* rest frame. The angular distributions [11] given below allow us to determine both 
the 1 I, I All I, I A_\_ I magnitudes and phases 6\\,5_\_ (up to a two-fold ambiguity [pX|, 12]). Thus 



1 (^3p q / X 

(6) 



where 



/i = 2cos^6'x*(l -sin^6'trCos^$tr) , /2 = sin^ ^'i^'. (1 - sin^ 6'tr sin^ ^> 



tr) , 



fs = sin 9k* sin 6tr , /4 = ± sin Ok* sin 2dtr sin ^ 



tr , 



1 2 1 

/s = y= sin 29k* sin Otr sin 2$tr , /e = ^—7= sin 29k* sin 26'tr cos <I>tr • 

\/2 \J2 
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The plus sign in /4 g refers to the B mesons which are bq (q = u, d) bound states, and the minus 
sign to B mesons (bq bound states). In the following, for convenience, the amplitudes are implicitly 
written for the B mesons, since we are dealing with the b quark and not the antiquark b. 

2-2 Effective Hamiltonian 

The basis for nonleptonic weak decays of hadrons is the operator product expansion, and the effective 
Hamiltonian relevant to B— > J/^' + K* may be written as 

1^ 6 

Heff = ^(^cfeK: [Cl{^^)Ol{^^) + C2{^i)02{^A] - Vt,v:,Jy,{^i)o,{^,)) , (7) 

where the Wilson coefficients are evaluated at next-to-leading order and at the renormalization 

scale jjL. We neglect the electroweak penguin operators O7...10 since their corresponding coefficients 
C7 • • • Cio being proportional to = 1/137 are numerically negligible compared to the dominant 
Ci and C2 associated to the tree diagrams, and C3 • • • Cg associated to the gluonic penguin loop 
diagrams. We have: 

V-A current x V-A current 

Oi = {SaCa)v-A{ci3hp)v-A , ©2 = {SaCi3)v-A{ci3ba)v-A , (8) 

QCD-Penguins 

O3 = {Saba)v-A^iqfiq/3)v-A , O5 = {Saba)v-A^{qf3<ll3)v+A , (9) 

1 g 



04 = {Sabf3)v-A^iq0a)v-A , Oq = {Sabf})v-A^{q/3qa)v+A , (10) 
9 9 

where a, (5 are quark color indices and in (9)-(10), the sum is done over q = u,d, s, c, b quarks. 
Also {qiq2)v-A{q3q4)vTA denotes q^jp{l - 75)^2 ^7^(1 T 75)94- The coefficients Ci(/u) are given 
in Table XXII of [|l^] at next-to-leading order in the naive dimension regularization (NDR) and in the 
't Hooft-Veltman (HV) 75 renormalization schema : 

Ci = 1.082 (1.105) , C2 = -0.185 (-0.228) , C3 = 0.014 (0.013) , 

C4 = -0.035 (-0.029) , C5 = 0.009 (0.009) , Cq = -0.041 (-0.033) , 

where the first numbers refer to the NDR scheme and those in the parentheses to the HV scheme, both 
evaluated at the scale = Wi^{nH)) = 4.4GeV. The dependence of Ci{p) on the scale ^ as well as 
on the regularization-schema must be cancelled in principle by the matrix-elements ( K*^'|Oj(^)|B ) 
since physical amplitudes ~ Ct(/i) x ( K*^'|C'j(/x)|B ) must be scale and regularization-scheme inde- 
pendent. In the "naive" factorization approach, ( K*^'|C'j(/i)|B ) is a product of decay constants and 
form factors, both are real, moreover they are scale and regularization scheme independent, hence the 
amplitudes via Ci{ii) suffer from these dependences and turn out to be ambigous. After a tentative 



approache||14|], the QCD methods finally solve this problem as we will see. 

Using the unitarity condition VubV*s + VchV*s + VtbV^g = and neglecting VubV*^ ~ 10~^, we 
note that (7) has only the unique VchV*^ factor. This is crucial to ensure that no matter the penguin 
( Ozfi ) and tree ( Oi^2 ) matrix-elements are, we always have | A4/A^ | = 1 (to a very good precision 
VubVus ~ 10"^), where M and JA are respectively the B^ J/^' + K(K*) and B J/^ + K(K*) 
amplitudes. This condition |A4/A^| = 1 allows us to extract from experimental data|15, [T^ the 



CP asymmetry j3 angle without any theoretical uncertainties[17] through the "gold-plated" modes 
B^ J/^- + K(K*). 
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2.3 QCD-improved factorization approach 

To proceed further, we use the QCD-improved factorization approach according to which, in the 
infinite b quark mass limit and for some classes of two-body hadronic Mi + M2 decays, the mass- 
singularities (infrared divergences) factorize, so that the amplitudes may be written as convolutions 
of universal quantities (the light-cone distribution amplitudes (LCDA) of the mesons B, Mi, M2 
with their associated semi-leptonic form factors F^^^ {rn^)) and QCD perturbatively calculable hard- 
scattering kernels {u),T^^ r], u) which are process-dependent. 

Here Mi denotes the recoiled meson which can be light (vr, p K, K* • • • ) or heavy (charm D, D* 
• • • ) but the emitted M2 can only be a light meson[^ or a QQ quarkonium (heavy or light) and not a 
heavy meson Uke D, D* • • • . The main problem with the latters is that they represent an extended soft 
hadronic object. In our case Mi is the K* and M2 is the J/'f. Schematically we write 

(k*^| |b) = F^K*(m|) y'duTi(u)cI>^(u) + j' d^dr]dnTf{^,r,,u)<^B{0 ^k-(??)$*(u) , 

(11) 



B 



(r 



1 + r„< + 0{AQCD/mb) 



(12) 



where Oi are products of different currents in (8)-(10). 

These equations imply that in the —>■ 00 limit, "naive" factorization (corresponding to the first 
term 1 in (12)) is recovered in the absence of QCD for which T\u) is independent of u and T^^ = 
in (11). Since the b quark mass is large but finite, power corrections 0{Kqcd /frib) could be 
significant, specially in some particular cases (for instance in B^ tt + K for which the scale is not 
Aqcd but chirally-enhanced like rn^/{'ms + md)). However it must be noted that in the QCD- 
improved approach[||], these power corrections (generally associated with nonleading higher-twist 
LCDA) cannot be reliably computed since in many (but not all) cases, mass singularities again show 
up thus do not factorize. Keeping this fact in mind, we nevertheless calculate to order the QCD 
corrections to the naive factorization B^ J/^* + K* amplitude|^, the simpler case B^ J/^' + K of 
B— >V-i-P has been previously studied[|T^, |T^ within the same theoretical framework. 

We also mention another approach!^ called PQCD (perturbative QCD) according to which the double 
logarithms Sudakov suppression effects could regulate the mass singularities, hence power corrections 
and form factors may be perturbatively calculable but questioned in[^]. 

The symbolically written form factors F^^*(m|,) in (11) and the more explicit ^ K* |J7^| B ^ ones 
in (12) are in fact defined according to 



K*(ei,p) |57;.(1-75)MB(P) 



(r 



B(P)) 



Va - A 



where [C 



V^ = ie^,„^ er(p)P>' 



a « 2V(q2) 



M + niK* 



M + niK* 



m 



9m 



+2mK.Ao('z')^9^ , (13) 

q^ 

and q = P — p is the four momentum of the emitted J/^' meson. In the above equation, terms 
proportional to vanish when multiplied to the J/^* polarization vector ef^iq)- With 

( I |o ) = ( ^(e2,q) |C7^ c |0 ) = hm^ e*/{q) 
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where w 405 MeV is the J/^* decay constant extracted from the J/^* ^ ^+ + £ rate, the first 
term (^K* | J"^] B ^ ( ^ | ^'''lO ) on the right hand side of (12) is equal to 

L M + TTiK* M + niK* 

(14) 

When we compare (1) with (14) using (2), (7) and (12), then in terms of an overall common factor 

K = ^VaV:,Um^{M + m*^)Ai{ml) , (15) 

and the dimensionless constants a, b, c given in[Q] 

^ - mj. - ml 2M^Kl 2MK, 

" 2m^m^ ' m^m*(M + m^)2 ' ^ (M + m^)2' *• 

together with the form factor ratios x, y defined by[^] 

A,{ml) ' y A,{ml) ' 

we obtain the helicity amplitudes Ha : 

Hq = K{a — bx)ao , }i±i = k {±cy — l)a± . (18) 

Numerically we have a = 3.165,6 = 1.308, c = 0.436 and the dimensionful factor k (in KeV) is 
equal to 2.48 ^i(m|,). 

It remains the most involved computations of the polarization-dependent coefficients oa in (18) which 
are the nonfactorizable correction terms derived from the QCD-improved factorization frame- 
work to which the next section will be devoted. In the "naive" factorization method previously 
considered[^] without the QCD-improved approach, the oa was the process-independent BSW[22] 



coefficient 02 = C2 + Ci/3, and the penguin contributions Cs^e are neglected. As is well known, 
this "naive" factorization method suffers from a serious problem related to scale /i and regularization- 
scheme dependences of the Wilson coefficients Cj(^), while the amplitudes Ha are not. The inclusion 
of the corrections, as will be shown in (32)-(35), cures this problem. 

For B J/^' + K* decay, we remark that the original V-A left-handed property of the current 
57^(1 — 75)5 in Oi is reflected by the expression (18) of the helicity amplitudes H±i ~ ±cV{q^) — 
Ai{q^). This may be understood as the following: The —1/2 helicity left-handed s quark, emitted 
through the V-A current, picks up the spectator antiquark q (which has both ±1/2 helicities) to form 
the K meson. Thus the latter can only have —1 or helicity and not +1, since the s quark would 
maintain its —1/2 helicity unless final state interactions (FSI) are strong enough to flip it into a +1/2 
state. Therefore |H_i| ~ \cy + 1| would largely dominate |H+i| ~ \cy — 1| unless very strong FSI 
would reverse the order by making |a+| 3> |a_|. As we will see in section 2-5, this possibility is not 
supported by our calculations within the QCD-improved factorization approach. However the answer 
as always must come from the experimental side : whether or not |H_i| » |H_|_i| can only be settled 
by future measurements of the muon polarization inB— >J/^ + K follows by J/^' — > /i"*" + /x". 

The decay rate is obtained using (4) and (18) from which we get the normalized longitudinal ^0 and 
transverse Ay, Aj^ fractions measured by Argus|§], Cleo[^], CDF|0] and BaBar[^] collaborations 

\A^\2 j-0\\-bx? ^ ^^^^ 
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^ |2 ^ I Q+ + g- P + cVl a+ - g- p - 2cy (|a+p - |a-p) 
III 2S 



2 ^ I a+ + g- Pc^ + I a+ - g- P - 2cy (|a+p - la-p) 
±1 2S , I J 



where 



S = I So P(a - bxf + (I S+ |2 + I 5_ |2) (1 + c2 y2) _ 2cy (| S+ 



The phases 6o, , 6± of the Aq, Ay , A± are given by those of the ax since k, a, b, c, x, y are real. 

Provided that | a- P > | g+ P which is true from the first line of (40), and for nonzero Ai{q'^), a 
remarkable upper bound for the longitudinal fraction {Aq p can be derived using (19), no matter how 
are the finite form factors Ai{q'^), ^2(9^) < {a/b)Ai{q'^) = 2A2Ai{q'^) and V{q'^): 

^0 < — , P = p^-fo ' (22) 

the derivation of this upper bound can be easily obtained by considering the lower bound of the inverse 
1/1^0 p. Of course when all of the three ax are real and identical, we recover our old result[^] of the 
naive factorization method, as it should be: 

I 4 |2 - - ^^)^ < _noo /o^x 

[a — bxY + 2(1 + c^y^j + 2 



A,, 1 2 _ nA\ 

" ' (a - 6x)2 + 2(1 + c2?/2) ' 



A |2 ^ ^"^^^^ /o^N 

^ ' (a - 6x)2 + 2(1 + c2y2) ' ^^-"^ 



The latest experimental data for Aq, A\\ , Aj^ are|Q ^ 



^0 P = 0.597 ± 0.028 ± 0.008 , | A_l |^ = 0.160 ± 0.032 ± 0.036 



I ^11 P = 1 - I ^0 P - I ^± P = 0.243 ± 0.034 ± 0.033 . (26) 
From the experimental value of | P = 0.6, we derive a constraint on the ratio p using (22), 

' 5+ |2 + I a__ |2 2 2 



go p 



< -a^ = 6.6 



which is however far from being saturated (section 2-5 below). 

For the phases 5i, since measurements of the interference terms in the angular distributions are limited 
to Re( A\\Al), Im( A^A^) and Im(A_L^|), there exists a two-fold ambiguity[|l|, 



<-> (5_L ^ vr — (5_L, — 5\\ vr — {5±_ — <5||) • 
The phases quoted in radians are[Q, |l^] 

6±_ = arg( A±Al) = -0.17 ± 0.16 ± 0.06 (-2.97 ± 0.16 ± 0.07) 
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5j| = arg( A\\A*q) = 2.50 ± 0.20 ± 0.07 (-2.50 ± 0.20 ± 0.07) , (27) 

the numbers in parentheses correspond to the second solution due to the mentioned ambiguity. Con- 
sequently, from (26)-(27), one deducesljl^] either |H+i/H_i| = 0.26 ± 0.14 or |H_i/H+i| = 
0.26 ± 0.14, this ambiguity can be solved in the future by the J/^' + i~ lepton polarization 

measurements. 

The most important information we can draw from the measured nonzero phase of is that nonfac- 
torizable corrections to the "naive" factorization method must be taken into account. 

2-4 Nonfactorizable Corrections 

In the QCD-improved factorization approach, the light-cone distribution amplitudes play a central 
role. For vector mesons, the LCDA are given by|^, |4|] 

{V{p,e) My)Qj{x) |0 ) = \l'^V ei(''P->-+(^-^)P-) (fyMy [^|$|| (r?)+ /^(^K??) 



fvMv 



(28) 



where i,j denote the Dirac spinor indices, z = y — x, ey* (e^) are the longitudinal (transverse) 

polarizations of the vector meson and fv and Fy are respectively its vector and tensor decay constants 
defined by 

{Vip,e)\q{x)-f^q{0)\0) = fvMv^p^. [ dr? e'^P-^.||(7y) , 

p-x Jo 



(29) 



Contracting the above equation ( V{p, e) \q{x)a^,yq{0) \ ) with p'^ and applying the equation of mo- 
tion together with the definition ( V{p, e) |g(0)7^g(0) | ) = fyMye^, a relation is obtained between 



the fv and Fl decay constants [|18|, |19|]: 



F^Mv = 2fvm^ 



fvMv 



m, 



-V 



Mv J 



2r/' 



(30) 



where in the last step the on-shell relation rf = has been applied! 18]. 

Finally are the leading twist-2 LCDA amplitudes while the nonleading vector-like 

twist-3 LCDA is denoted by (^]^(?/). The axial-like twist-3 LCDA 'p\{'q) contribution is negligible 
since it is proportional to fyMy — 2Fymq = My — Am"^ 0. We thus introduce expressions 
$y^*(r/), ^Yiv)' fT^iv) and ipY'^iv) for the vector meson K*. 
For the pseudoscalar B meson, its wave function in the heavy b quark limit is 



( B(P)|6i(a;)g,(0)|0> U+=,^=o 



iff 



(f + M)75$b(0 + 



(31) 



where • • • denote terms that do not contribute to the decay amplitude calculated later in (39)-(40). 

For the vector meson J/"^ treated as a heavy charmonium, to the leading order in l/rric, its wave 
function has a similar expression: 



( J/^ip,e) \ci{x)cj {0)\0) 



ki 
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We note that the use of the hght-cone wave function for the heavy J/^ is problematic, higher twist 
effects have to be included, however they may not converge fast enough. We adopt in the following 
<I>II (u) as the distribution amplitude (DA) of the nonlocal vector current of J/^, thus we treat the J/^f 
wave function on the same footing as the B meson. Comparing the above equation with (28), we see 
that at the leading order in l/rric one has for the heavy charmonium J/"^ : 

#11 (u) = ^±{u) = and = fq, . 

Equipped with these ingredients, we are ready to compute the nonfactorizable correction terms. 

Loop integrations of the four vertex corrections diagrams (Fig.6 in[||] which gives F^) and the two 
spectator diagrams (Fig. 8 in[||] which gives F^j) are not detailed here. Only we give the results in the 
NDR scheme: 

ax = (a^ + a^ + a^) , (32) 



«3 



«5 



C3 + ^ + 



4tt 


2iV2 


as 




47r 


2iV2 


as 





rub 



Nr 47r 2iV2 



-6 + 121n — + + F/^ , 



(33) 



(34) 



(35) 



with A'^c = 3 is the color number. For completeness, the constants —18, —18, —6 in (33)-(35) be- 
come respectively —14, —14, —18 in the HV scheme. These constants and the ln{mi,/n) term inside 
the parentheses of (33)-(35) reflect the scale and regularization-scheme dependences, they are can- 
celled by those of the Wilson coefficients Ci(/i), therefore the final expressions of ax are // and 
regularization-scheme independent. The F^ and F^j are calculated to be: 



Ff 



du 



$||(u)/C(r,u) +$x(u)W(r,u) 



where 



^ on N -1 3(l-2u)lnu 2r(l-u) 
/C(r, n) = 3 [ ln(l - r) - ivr] + -^-^ '- + ^ ^ 



+ 



l-u 



u 



(1-to)2 [l-r(l-u)]^ 



r nln(rn) + 



1 — u 1 — ru 

r'^u^l ln(l — r) — ivr] 



and 



H{r, u) = 8r 
with r = m\, /M^. Also 







1 — r(l — u) 1 — ru 
Jo 

r li 

^ (l-u) 



[l-r(l-u)]2 ' 

ln(l - r) -ivr-i 

m ru -, ^ > 

^ ' 1-r l-u J 



(36) 



(37) 



, ^ , . ( Inu rln(ru) ln(l — r) — ivr "i 

+8/ du$^uu2 + - ^^^-r / / 

L l-u 1 — r 1 — u 1 — r 1 — u J 
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»1 



J ^ , J (3-2u)lnu 2ru / 3 l-r(l+u) \ ^ , , 

du^ll -I + T+ "1 + 71 A ^ ruin ru 

" 1 — u 1 — r(l — u) V 1 — ru [1 — r(l — ujj^/ 



+ ( 3 - r(6 - 5n) + r^(3 - Stt + W)) .^"^^ — ^ I . (38) 



[l-r(l-u)P 



Similarly to the J /'^ +K case found in[|l8|, [19|], the infrared divergences in are mutually 



cancelled among the four vertex diagrams, this cancellation is essentially due to the symmetric u ^ 
1 — u of the kernels /C(r, u), Tl{r, u) and the wave functions. We note that in (38), the J/^ nonleading 
ip\{u) wave function also contributes to the transverse Ff on the same footing as the leading ^_l{u). 

For the F/y of the spectator-quark effect, following [20] we get 



Nc m^{M + mK*)Ai{m^){a -bx) Jo ^ Jo V Jo u 



N, M^M + mK.)^i("^|)(l T cy) Jq " ^ 



We first emphasize that the infrared-finite longitudinal Ff^ in (39) is not 1/M power-suppressed 
contrary to its appearance, since the B meson wave function $b(0 is appreciable only for ^ of the 
order Kqcd/M, hence the integral J d^<I>B(0/'^ ~ M/Aqcd compensates the fK*/{M + mx*) in 
(39). 

As for the transverse Ffj parts given in (40), they are infrared divergent although \/M and 
power-suppressed. Indeed, the first term in F^} ~ J (77) / 77^ unexpectedly diverges even with the 
K* leading twist-2 LCDA <I>^*(r/), the remaining divergent terms come from the twist-3 LCDA of 
both the K* and J/^' mesons. We have neglected the r dependences in to simplify the computa- 
tions of complicated loop integrals, since infrared divergences occur no matter the r dependences are 
kept or not. 

In the numerical applications, we use for the leading twist-2 LCDA their asymptotic form = 
^±(3^) = 6x(l — x) for both J/^* and K* mesons, although the former is treated as heavy. The twist-3 
LCDA of the K* vector meson are taken to be ip\{x) = 3[1 + (2x - l)^]/4 and ip\{x) = 6x(l - x). 
For the B meson wave function, we take ^>b(0 = - 0^ exp[-S,'^A'P /2u;'^] with to = 0.25 

GeV, and Nb is the normalization factor such that jJ^d^^BCO = 1- With this <^b(0' we get 
Jq ^b{0/^ = M/(0.3GeV), in agreement with our guess / dS,^B{0/^ ~ M/Aqcd mentioned 
above. 

Remarks 

From (36)-(40) we draw some unexpected features of the nonfactorizable terms in B^ J/^' + K* 
which are dictinctive from B^ J/*? + K : 

1- The chirally-enhanced factor m^/ {ms+mii) inherent to the pseudoscalar K meson inB^ J/^f+K 
is absent in B— > J/^' + K* with the vector K* meson. Therefore, neither Ff nor F^^ are chirally- 
enhanced in both the longitudinal and transverse parts. 
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2- Both the l/rric leading and nonleading wave functions of the charmonium J/^' contribute to the 
transverse and Pf^, moreover Pf is infrared-finite with the nonleading t^\{u) as shown in the 
first line of (38). The longitudinal Pf and Pfj are also infrared-finite and not power-suppressed. 

3- At the first order Aqcd/M level, P^j vanishes, only survives the infrared-divergent Pfj. Unex- 
pectedly, the infrared divergence of Pfj already comes from the leading twist-2 LCDA of the K*, its 
twist-3 is unnecessary to render Pfj divergent. At the second order Ag^^^j/M^, the twist-3 ip^* {t]) 
and (P_l{u) of both K* and J/^ respectively make P^j infrared divergent. 

4- Fortunately, the longitudinal part clq as given by (36) and (39) is infrared convergent, therefore it 
can be unambigouly used in the following section to check whether or not agreement exists between 
theoretical calculations and experimental data. 

2-5 Numerical Results 

In (39) we need the decay constant /k* of the vector meson K*, for that we may use the r lepton 
decay r — > i/,- + K* width given by 



to extract /k* ~ 210 MeV, a value consistent with mK*/K* = i^^pfp of the SU(3) flavor symmetry 
which relates /k* to the decay constant fp w 198 MeV of the charged p(770) vector meson. We also 
take /vi, w 405 MeV and /b 180 MeV. 

Numerical values of the infrared-finite quantities Pj in (36)-(37) and Pfj in (39) are 

i70 = _o.82 - i 6.61 , = , , 2^:7 ■ . . (41) 

Ai(m|)(a - bx) 

These numerical values may vary within 20 per-cent when different LCDA are used instead of their 
asymptotic forms. Since both Pj and Pfj are finite, we first concentrate on the longitudinal part 

where Br(B J/^' + K*) is the branching ratio = (1.35 ± 0.18) x 10"^ for the neutral and (1.47± 
0.27) X lO^'^ for the charged B meson and tb is their respective hfetime [1.56(1.65) ±0.04] x 10^^^ s. 
We obtain on average rL(B J/^- + K*) = (3.37 ± 0.4) x lO^^^ GeV, using |^oP = 0.6. When 
this experimental value is compared with the theoretical expression 

Fl(B ^ J/^ + K*) = ^^^\a - bxflZof , 

we find that the product (a — 6x)^i(m|,)|ao| is constraint to equal 0.156 ± 0.02, thus 

(a - 6x)yli(m|)|So| = 0.156 ± 0.02 . (42) 
Our formulae (33)-(39) with as(mh) = 0.23 give |ao| ~ 0.14. 

This value of |ao| in turn can easily satisfy the constraint (42), and we get a domain for x, Ai(m|,) 
plotted by a hyperbolic curve and translated into the following numerical values: 

< X < 1.1 and 0.35 < Ai(m|) < 0.60 . (43) 

The smallest x is associated with the smallest Ai(m|,), the latter increases with increasing x. For the 
transverse a± which cannot be reliably calculable because of their infrared divergences, we reverse 
the naive factorization procedure previously proposed[^] in which (23) was used to determine x,y. 
Now we fix X 1.1 and |ao| = 0.14 then using the theoretical expressions (19)-(21) matched with 
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the experimental data (26)-(27), we determine in turn a± and y. The resulting contour solutions 
|a+| 0.095 ± 0.02, |a-_| 0.125 it 0.02 confirm the polarization-dependence of ax- Also we get 
y « 1.75. 

We remark that our favoured values x = A2{m'^)/Ai{m'^) < 1.1 and y = V {m^) / Ai{m'^) 1.75 



,27 



. It 



are generally satisfied by some models of form factors studied in the literature[^ 
is amusing to note however that x < 1.1 is at variance with the B-K* form-factor ratio derived below 
from equation of motion for on-shell massless strange quark. Indeed using 



P'^s-f^j^b = russ-f^b = , 
then we obtain a relation between the two form factors Ai(g^) and A2{q'^). Assuming 

(K*(ei,p) I 57^756 I B(P) ) = = , 



(44) 



which gives 



we get 



(M + mf{M'^ -m^ - q^^ - 
\{M\m\q^) 



2m Ao{q'^ 



M + m^i(g2) 



(45) 



From (45), we recover the well-known relation at g2 = 



M + m 



MO) 



M -m 



^2(0) = Ao{0) . 



Neglecting m"^ with respect to M"^ and q"^ = m^, then we get from (45) 

A2{ml) _ M2 



X = 



Ai{ml) M2-m| 



1.52 , 



this too large value indicates that = in (44) may not be a good approximation. 

Although the relation (45) is derived here with the assumption nis = 0, we remark nevertheless that it 
bears some similarity with the one [28] derived in the very different context of the large recoil energy 
g2 hmit for which m should be neglected: 



A2iq' 



(M + m) 



^i(^2) M2-m2-g2L M2 



1 



2Mm 



M2 



M2 



^2 ■ 



This is only this large q"^ recoil energy limit that can justify the above result of [28] and not the rris 
assumed here in (44). 



3 Conclusion 

We have examined within the QCD-improved factorization approach different aspects of the color- 
suppressed B decay into two vector mesons B^ J/^ + K* for which important experimental results 
are recently obtained[^ |5|]. The nonzero phases of the helicity amplitudes measured for the first time 
by[§, ||] indicate that nonfactorizable terms must be taken into account. We emphasize that the phases 
can only be determined in B^ V -1- V decay, hence its superiority over the B^ V -1- P and B^ P -1- P 
in this aspect. Our calculations (36)-(39) give nonzero imaginary part to the process-dependent and 
polarization-dependent coefficients ax{"^K*) that substitute the conventional process-independent 
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and polarization-independent BSW 02 coefficient. The phases obtained for our ax{'^K*) are consis- 
tent with experiments. 

We find that the longitudinal ao(^'K*) is free of infrared divergence and |ao(^K*)| ~ 0.14 is about 
twice the BSW 02 ~ 0.074, thus corrections - mainly due to the spectator-quark effect Fjj in (38) 
- are large but under control. This |ao(^'K*)| is also different from the one in B^J/^' -i-K case for 
which experimental data indicate that |a2(^'K)| 0.25, thus confirming their process-dependence. 

On the other hand, our calculations show that the transverse part a±{^K*) is infrared divergent 
(although power-suppressed), this infrared divergence may be handled by a cutoff procedure. From 
remarks in 2.4, we note an important difference between B^ J/'f + K and B^ J/'f + K* in their 
nonfactorizable terms. In the former case, the discrepancy by a factor of three between experimental 



data (lai(^'K)l 0.25) and the theoretical estimates [Il8|,|l9|] (|a^(^'K)| 0.08) using twist-2 LCDA 
may be cured [|T8l] by introducing the chirally-enhanced twist-3 LCDA of the K meson which gives 
a formally infrared divergent a2(^K)- This divergence is parameterized[Q, [T^] by a random number 
X just to render 02 (^K) finite and large enough to fit data. This input cannot be evoked here for the 
B^ J/^f + K* case, since the longitudinal ao(^'K*) (the equivalence of 02(^'K)) is finite. Moreover, 
if one postulates that the input X is used to cure the infrared divergence of a-|-(^'K*) - similarly to 
the B ^ J/^f + K case where X cures the 02 ( ^K) - then the experimental data on | ^ 1 1 1 , | A± \ < \Aq\ 
imply that the input X is constraint to make the divergent a±(^'K*) be smaller than the convergent 
ao(^'K*), which is somewhat disturbing. 

Therefore we are inclined to believe that in our case of B— + K*, the procedure used to handle the 
infrared divergence via X may be inadequate for the treatment of the discrepancy (if any) between 
experimental data and theoretical estimates using QCD-improved factorization approach. We may 
seek the remedy outside the ax, probably in the form factor Ai (m|,) and in the ratios x = A2/A1 and 
y = V I A\, since the overall factor A\(m^) in (15) is central to the absolute strength of the decay 
rate B— > J/^f + K*, as well the ratio x is central to the longitudinal part and y to the transverse 
fraction A^. To deal with the transverse part a±(^'K*), we adopt a pragmatic method by fixing x 
and |ao(^'K*)| previously obtained from \Aq\, then using data on |j4||| and \ Ax\ together with their 
theoretical expressions (19)-(21), we determine in turn |a+(^'K*)| « 0.095 ± 0.02, |a_(^'K*)| k. 
0.125 lb 0.02. Moreover the ratio y = V{'m\,)/Ai is also bounded around 1.75. The constraints 
X < 1.1 and y k, 1.75 have implications on models of B-K* form factors commonly used in the 
literature. 

In summary, our results show that the spectator effects and final state interactions reflected by F'^^ 
play an important role in our quantitative understanding of the color-suppressed B— > J/^' -i-K* decay 
for the dominant longitudinal mode. However the power AgcD/m-b corrections within the QCD- 
improved factorization approach has to be better understood. 



Note added: After the first version of this paper is circulated, we received a preprint |hep-ph/01 11094 
by H.Y. Cheng, Y.Y. Keum and K.C. Yang. The second version takes into account their greatly 
appreciated remarks. We thank the referee of JHEP for her (his) critical remarks, helping us to 
improve the paper in its third version. 
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